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Abstract 

We have proposed a model for relativistic compact star with anisotropy and analyt¬ 
ically obtained exact spherically symmetric solutions describing the interior of the 
dense star admitting non-static conformal symmetry. Several features of the solu¬ 
tions including drawbacks of the model have been explored and discussed. For this 
purpose we have provided the energy conditions, TOV-equations and other phys¬ 
ical requirements and thus thoroughly investigated stability, mass-radius relation 
and surface redshift of the model. It is observed that most of the features are well 
matched with the compact stars, like quark/strange stars. 
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1 Introduction 


Since the striking idea of white dwarf hj Chandrasekhar [T] the study of general 
relativistic compact objects received a tremendous thrust to carry out research 
in the held of ultra-dense objects. White dwarfs are composed of one of the 
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densest forms of matter known, surpassed only by other compact stars such 
as neutron stars, quark stars/strange stars, boson stars, gravastars etc. 

In the compact stars the matter is found to be in stable ground state where 
the quarks are conhned inside the hadrons. If it is composed of the de-con£ned 
quarks then also a stable ground state of matter, known as ‘strange matter’, 
is achievable which provides a ‘strange star ’ iZEEE] . The main theoretical 
motivation for postulating the existence of strange stars was to explain the 
exotic phenomena of gamma ray bursts and soft gamma ray repeaters PI7] . On 
the observational side, the Rossi X-ray Timing Explorer has now conhrmed 
that S'y4XJ1808.4 — 3658 is one of the candidates for a strange star [8]. 

As a start of the astrophysical objects of compact nature of the above kind 
the theoretical investigation has been done by several researchers by using 
both analytical and numerical methods. However, the density distribution in¬ 
side the compact stars need not be isotropic and homogeneous as proposed 
in the TOV model. In the early seventies Ruderman |9] has investigated the 
stellar models and argued that the nuclear matter may have anisotropic fea¬ 
tures at least in certain very high density ranges (> lO^^gm/c.c.), where 

the nuclear interaction must be treated relativistically. Later on Bowers and 
Liang [10] gave emphasis on the importance of locally anisotropic equations 
of state for relativistic fluid spheres. They showed that anisotropy might have 
non-negligible effects on such parameters like maximum equilibrium mass and 
surface redshift. Anisotropic matter distribution have been considered recently 
by several investigators as key feature in the conhguration of compact stars 
(with or without cosmological constant) to model the objects physically more 

As one of the physical characteristics Mak and Harko (T^j have calculated the 
mass-radius ratio or compactness factor for compact relativistic star however 
in the presence of cosmological constant. To study mass and radii of neutron 
star Egeland [20| incorporated the existence of cosmological constant propor¬ 
tionality depending on the density of vacuum. Kalam et al. [18] by using their 
model calculated mass-radius ratio for Strange Star Her X — 1 oi their model 
and found it as Mptr/R = 0.336 which satishes Buchdahl’s limit and corre- 
sponds to the same mass-radius ratio for the observed Strange Star Her X — 1. 
This physical aspect is very important because of the fact it helps to deter¬ 
mine nature of the compact stars, its size, shape, matter contain and several 
other features. 

Now to get more information, several studies have been done on charged or 
neutral fluid spheres with a spacetime geometry that admits a conformal sym¬ 
metry, in the static as well as generalized non-static cases. In this line of 
conformal symmetry we note that there are lots of works available in the lit¬ 
erature PT]|22II23II24II25| . However, we investigate in the present work a new 
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anisotropic star admitting non-static conformal symmetry (mathematical for¬ 
malism is provided in detail in the next Sec. 2). 

Under this background and motivation we investigation in the present pa¬ 
per a model for relativistic compact star with anisotropy and hnd out exact 
spherically symmetric solutions which describe the interior of the dense star 
admitting non-static conformal symmetry. The scheme of this study is as fol¬ 
lows: We provide mathematical formalism for non-static conformal motion in 
Sec. 2 whereas the Einstein held equations for anisotropic stellar source are 
given in Sec. 3. The solutions and general features of the held equations under 
the Dev-Gleiser energy-density prohle have been found out in Sec. 4 along 
with a special discussion on two prescriptions, in section 5, viz. {5.1) the Mis- 
ner prescription, and {5.2) the Dev-Gleiser prescription are done to show that 
mass to size ratio clearly indicates that the model represents a quark/strange 
compact star. In Sec. 6 we have made some concluding remarks. 


2 NON-STATIC CONFORMAL SYMMETRY 


Now, we consider the interior of a star under conformal motion through non¬ 
static Conformal Killing Vector as 


^i9ij 9ij-,k^ T 9kj^.^i T 9ik^.j '4^9iji 


( 1 ) 


where L represents the Lie derivative operator, ^ is the four vector along 
which the derivative is taken, -0 is the conformal factor and gij are the metric 
potentials |26]. It is to be noted that for 0 = 0, the equation above yields the 
killing vector, whereas for 0 = const, corresponds to the homothetic vector. 
So, in general for 0 = 'ijj{x,t) we obtain conformal vectors. In this manner, 
one can perform a more general study of the spacetime geometry by using 
Conformal Killing Vector. 

The static spherically symmetric spacetime is given by the line element (with 
G = 1 = c in geometrized units) [3T1I32] 


-I- -|- r‘^{d9‘^ + sin^Odcj)^), 


( 2 ) 


where i>{r) and A(r) are the metric potentials and functions of radial coordi¬ 
nate r only. 

The proposed charged fluid spacetime is mapped conformally onto itself along 
the direction Here following Herrera et ah [331134] we assume ^ as non-static 
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but ^|J to be static as follows: 


^ = a{t,r)dt + /3{t,r)dr, 

( 3 ) 

Ip = ip{r). 


( 4 ) 

The above set of eouations f 11-141 eive the followine: exoressions [261127112811291115 
for a, f3, Ip and A: 

a = A + i^kt, 


( 5 ) 

pj = -Bre~^, 

2 


( 6 ) 

Ip = Re“2, 


( 7 ) 

u /^2 2 

e = C r exp 

A 

1 /* 62 

-2kB-^ / —dr 
J r 

( 8 ) 


where k, A, B and C are arbitrary constants. According to Maartens and 
Maharaj izg one can set A = 0 and B = 1 so that by rescalling we can get 


a = -kt, 

13 = -re 2 , 

_ A 

Ip = e 2 , 


( 9 ) 

( 10 ) 

( 11 ) 


A 



( 12 ) 


However, one can note that the above considerations, i.e. A = 0 and B = 1, do 
not loss any generality. This is because the rescalling ^ and pj in the following 
manner, p —)• B~^p and ip B~^ip, leaves Eq. (1) invariant. 


3 THE FIELD EQUATIONS FOR ANISOTROPIC STELLAR SOURCE 


In the present investigation we consider the most general energy-momentum 
tensor compatible with spherically symmetry in the following form: 

T!P = (P + Pr)u^U^ + PrgiP + {pt - Pr)V^Vu, (13) 
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with = —ri^r]fj_ = 1 , where p is the matter density, pr and pt are radial 
and the transverse pressnres of the flnid respectively. 


For the metric given in Eq. ([2]), the Einstein held eqnations are [35] 


e 


A 




e 


A 





Snp, 

(14) 

00 

(15) 



X'ly' — {u'— X') = STTpt. 

(16) 


Imposing the conformal motion, one can write the stress energy components 
in terms of conformal fnnction as follows: 


Snp = ^(1 — — 2ripip')^ 

(17) 

Snpr = — ‘^kip — 1), 

(18) 

Snpt = ^ {p) — kY + 2r'ip'ip' . 

(19) 


4 THE SOLUTIONS AND FEATURES OF THE FIELD EQUA¬ 
TIONS FOR DEV-GLEISER ENERGY-DENSITY PROFILE 


Now onr task is to hnd ont solntions for the above set of the modihed Einstein 
eqnation in terms of the conformal motion. However, to do so we wonld like 
to stndy the following case for specihc energy density as proposed by Dev and 
Gleiser for stars with the energy density in the form 

+ /Stt, ( 20 ) 


where a and b are constants. It is to be noted that the choice of the valnes 
for a and b is dictated by the physical conhgnration nnder consideration, e.g. 
a = 3/7 and 6 = 0, corresponds to a relativistic Fermi gas for nitradense 
cores of nentron stars [SB] whereas the valnes a = 3/7 and 6 7 ^ 0 represent 
a relativistic Fermi gas with core immersed in a constant density backgronnd 
[2T] . It can easily be observed from Eq. (|20|) that for large r the constant 
density term dominates (r^ 3> a/36) so that the corresponding star can be 
modeled with shell-type featnre snrronnding the core. 
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Due to use of the above energy density, Eq. ffT7|) will take the following form: 


+ 3 ^) = ^ ’ ( 21 ) 

Solving this Eq. (l2l]) . we get 

^jJ‘^ = 1 — a — br'^ + - = g{r) (say), (22) 

r 


where d is an integration constant. 


Substituting the value of -0 in Eqs. flT8|l and flT^ . expressions for radial and 
tangential pressures can be obtained as 


Pr = 


Pt = 


Svrr^ 

1 

Svrr^ . 


2 — 3a — 3br^ H- 2kJg(r 


1 — a — 3br^ + — 2k^Jg{r) 


(2.3) 

(24) 


y = 0,4' =0.15' W = 0.2’ 



Fig. 1. Variation of radial pressure, pr, with respect to the parameters a and b for 
a fixed radial distance r = 0.4 km 

In Figs. 2 - 3 we have shown the behaviour of radial and tangential pressures 
and whereas Fig. 1 represents the 3-dimensional presentation indicating the 
variation of radial pressure with respect to the parameters a, b for a hxed 
radius. All the curves are in accordance to the physical feasibility within the 
acceptable range. 

Let us now consider a simplest form of barotropic equation of state (EOS) 
as p = ijJip, where Wj are the EOS parameters along radial and transverse 
directions. Here at present as such we are not considering different forms of u 
with it’s possible space and time dependence as are available in the literature 
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'a = 0,85' 'h = 0,01’ 'k = 0.15"rf = 0,2’ 



Fig, 2, Variation of radial pressure with respect to radial distance 


'a = 0,95’ 'b = 0,01' 'k = 0.08’ 'd = 0.2' 

0.3 0,4 0.5 0.6 



Fig. 3. Variation of transverse pressure with respect to radial distance 

[37|38f39|l4n] . Thus, the EOS parameters u can straightly be written in the 
following forms: 



1 

(a + 36r^) 


2 — 3a — 36r^ H- 2k\J g{r) , 


(25) 



1 

(a + 3&r2) 


1 — a — 36r^ + — 2k\J g{r) 


(26) 


In Figs. 4 and 5 we plot variations of radial and transverse EOS parameter with 
respect to radial distance. In which radial EOS represents expected nature with 
a decreasing pattern. However, transverse EOS parameter indicates that it lies 
between (—|, 0). It is negative like exotic matter in natnre, however, later we 
will see that all energy conditions are satished for real matter situation. 
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'a = 0,85' 'h = 0,01’ 'k = 0.15"rf = 0,2’ 



Fig, 4, Variation of radial EOS parameter with respect to radial distance 


'0-0.95' 'i -0,01’ '* -0.08' 'rf-0,2’ 



Fig, 5, Variation of transverse EOS parameter with respect to radial distance 
4-1 Anisotropic behavior 


For the model under consideration the measure of anisotropy in pressures can 
be obtained, from Eqs, (23) and (24), as follows: 


A = {pt- Pr) 


1 

STrr^ 


; 2 o 3d 
h -|- 2ci — — 


(27) 


It can be seen that the ‘anisotropy’ will be directed outward when pt > Pr 
i,e, A > 0, and inward when pt < Pr i.e, A < 0, This feature is obvious from 
Fig, (5) of our model that a repulsive ‘anisotropic’ force (A < 0) allows the 
construction of a more massive stellar distribution |TT] , 
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'a = 0,95' 'h = 0,01’ 'k = 0.08' 'd = 0,2’ 



Fig, 6, Anisotropic behaviour at the stellar interior with respect to radial distance 
4-2 Energy conditions 

We now put and verify different energy conditions of the stellar configuration 
as follows: 

P > 0, P + Pr>0, P + Pt>0, P + Pr + 2pt>0, 


p > \pr\ and p > \pt\. 

It is seen that Null Energy Condition (NEC), Weak Energy Condition (WEC), 
Strong Energy Condition (SEC) and Dominant Energy Condition (DEC) i,e, 
all the energy conditions for our particular choices of the values of mass and 
radius are satished which can also be observed from Fig, 7, 


’fl = 0,85’ ’6=0,0r '/t = 0.15’ 'd = <d.T 



Fig, 7, Energy conditions as a function of r are plotted for the specified range 
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4^.3 TOV Equation 


To search equilibrium situation of this anisotropic star under different forces, 
we write the generalized Tolman-Oppenheimer-Volkoff (TOV) equation as [H] 


Mcip + Pr) dpr 2 

-5-^ + -(Pt - Pr) = 0, 

dr r 


(28) 


where Mq = Mg{t) is the effective gravitational mass inside a sphere of radius 
r which can be derived from the Tolmam-Whittaker formula as 


Mc^r) 


-L 2 — / 

—re'^u. 


(29) 


The above equation explains the equilibrium conditions of the fluid sphere due 
to the following hydrostatic, anisotropic and gravitational forces: 


Fh = 


dpr 

dr 


1 

Svr 


4 6a 9d k{2br‘^ + J) 


7"3 ^3 ^4 


9{r) 


Fa = -{Pt-Pr) = ^ 
r Sir 


4a 2k^ 6d 


rpO /^O 


(30) 

(31) 


V 


Fg = -^(P + Pr) 


1 

Stt 


2 , 12a 2 15d 10PVc/(r) l9akJg{r) isd? 6a^ 

{ -^ ^ ^- \\ ^ 


9{r)r 


15ad 


/yt V /yO ry ^ ryOny‘ 


2r^ 




The Eq. (28) can be rewritten in the form 

Fg + Fh + Fa = 0. 


(33) 


The prohles of different forces Fg, F^, Fa are shown in Fig. 8. The hgure 
exposes that our model of anisotropic star is in static equilibrium under the 
gravitational {Fg), hydrostatics {Fh) and anisotropic {Fa) forces. 
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'a = 0,85' 'b = 0.01’ 'k = 0.15’ 't/ = 0,2’ 



Fig, 8, The dark energy star is in static equilibrium under gravitational 
(Fg),hydrostatics (F/j) and anisotropy {Fa) forces, 

4-4 Stability issue 


The velocity of sound should follow the condition = dp/dp < 1 for a phys¬ 
ically realistic model [^2llT8] , We therefore calculate the radial and transverse 
sound speed for our anisotropic model as follows [H]: 


2 dpr 1 


4 — 6a -k — — 4:k\Jg{r) 


k{2br‘^ + 


9{r) 


2 dpt 1 
dp 2a 


2-2a + 2k^ - ikJgip) - ^ ’’ 


(34) 

(36) 


'fl = 0,85' F =0,01’ '/t =0.15’ 'rf = 0,2’ 



Fig. 9. Variation of radial sound velocity with respect to radial distance. 

Let us check whether the sound speeds he between 0 and 1. For this require¬ 
ment we plot the sound speeds in Figs. 8 and 9. It is observed that numerical 
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'a = 0,85' 'h = 0,01’ 'k = 0.08' 'd = 0,2’ 



Fig, 10, Variation of transverse sound velocity with respect to radial distance, 

'a = 0,85' ’/> = 0,01’ 'k = 0.08' 'd = 0,2’ 



Fig, 11, Variation of \v‘^J — with respect to radial distance, 

values of these parameters satisfy the inequalities 0 < \vlg\ < 1 every¬ 
where within the stellar object. However, for radial sound velocity there is a 
prominent change over after attaining certain higher value which gradually 
goes down. Since the transverse sound velocity is negative (due to negative 
transverse pressure), we use numerical value. As sound speeds he between 0 
and 1, we should have \v^g\ — v^g < 1 as evident from Fig, 10, 

Now, one can go through a technique for stability check of local anisotropic 
matter distribution as available in the literature |12] , This technique is known 
as the cracking concept of Herrera and states that the region for which radial 
speed of sound is greater than the transverse speed is a potentially stable 
region. Fig, 10 indicates that there is a change of sign for the term \v‘^g\ — v^g 
within the specific configuration and thus confirming that the model has a 
transition from unstable to stable configuration. The present stellar model 
gradually gets stability with the increase of the radius. However, in terms of 
the maximum allowable compactness (mass-radius ratio) for a fluid sphere as 
given by Buchdahl |1S], the stability issue will be discussed later on. 
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4-5 Surface redshift 


One can calculate the effective gravitational mass in terms of the energy den¬ 
sity p as 


R 

Meff = 4:71 j pr'^dr 
0 


-R 


a + bR^ 


(36) 


'q=0.85’ '6 = 0.01' 



Fig. 12. Mass vs radius relation is shown in the plot for the specified range 
Therefore, the compactness of the star is given by 


M. 


u = 


eff 


R 


= i [a + bR^ 


(37) 


The nature of variation of the above expression for compactness of the star 
can be seen in the Fig. 9 which is gradually increasing. We now define the 
surface redshift corresponding to the above compactness as 


l-b = [1 - 2 m] 



(a + bR^) 



(38) 


so that the surface redshift can be expressed as follows: 




1 - (a + bR^) " 


- 1 . 


(39) 


We plot surface redshift in Fig. 9. Likewise compactness factor u this has also 
a behaviour of gradual increase. This feature is also observed from the Table 
1. The maximum surface redshift for the present stellar configuration of radius 
0.62 km turns out to be Zg = 0.35280 (see Table 1). 
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'q=0.85’ '6 = 0.01' 



Fig. 13. Variation of redshift 2: with radial coordinate r 

Table 1 

Calculation of masses and hence compactness factors for different values of constant 
b with a = 3/7 and R = 0.62 km (with the conversion 1 km = 1.475 Mq) 


b 

0.0 

0.01 

0.02 

0.03 

Me// (in Mq) 

0.19596 

0.19952 

0.20335 

0.20717 

u 

0.21428 

0.21817 

0.22236 

0.22653 

Zs 

0.32279 

0.33257 

0.34257 

0.35280 


5 SPECIAL PHYSICAL ANALYSIS OF THE MODEL BASED 
ON MATTER CONTAINED 


If we look at Eq. (35) as well as Eq. (36) then it will at once reveal that 
there are lots of information inherently hidden inside these two equations. Let 
us therefore examine the present model for the specihed values of a and b, 
appearing in Eq. (35) and Eq. (36), as follows: 


5.1 The Misner prescription 


For the prescription of Misner [36] i.e. a = 3/7 and 6 = 0, we get the effective 
gravitational mass as M^ff = (3/14)i? which numerically comes out to be 
0.195 Mq so that the compactness factor becomes u = M^ff/R = 0.21428 as 
can be seen from Table 1. Here to estimate mass we have taken radius of the 
star R = 0.62 km by solving pr = 0 at r = R graphically from Fig. 11. This 
therefore represents a relativistic Fermi gas for ultradense cores of neutron 
stars as noticed in the Ref. [SB] . 
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'fl = 0,85' ’ft = 0,01"A=0.08’ 



Fig. 14. Radius of the star where pr = Q increases with the increase of the parameter 
d 

5.2 The Dev-Gleiser prescription 


In this prescription [2T], by snbstitnting the valnes a = 3/7 and 6 7^ 0 
one can get the expression for the effective gravitational mass as M^ff = 
(1/2)7? [(3/7) + bR"^]. For this case we provide here a data sheet by choosing 
different values for b (for Case A, 6 = 0 whereas for Case B, 6 > 0) in Table 1. 
This corresponds to a relativistic Fermi gas with core immersed in a constant 
density background [2T] . 


6 CONCLUDING REMARKS 


In the investigations we have observed some interesting points which are as 
follows: 

(1) Note that though transverse pressure is negative, however, all energy con¬ 
ditions are satished. This indicates that matter distribution of the star is real 
i.e. star is comprising of normal matter. 

(2) It is observed that the transverse sound velocity is negative (due to nega¬ 

tive transverse pressure) so that we use numerical value. As sound speeds he 
between 0 and 1, we should have < 1 (Fig. 10). 

(3) In connection to isotropic case and in the absence of the cosmological con¬ 
stant it has been shown for the surface redshift analysis that Zg <2 [lSll46ll47j . 
On the other hand, Bohmer and Harko m argued that for an anisotropic star 
in the presence of a cosmological constant the surface redshift must obey the 
general restriction Zg < 5, which is consistent with the bound Zg < 5.211 as 
obtained by Ivanov [18]. Therefore, for an anisotropic star without cosmolog- 
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ical constant the above valne Zg = 0.3 seems to be within an acceptable range 
(Fig. 12). 

(4) The radius of the star depends on the parameter d, as appears in Eq. (23), 
with an increasing mode (Fig. 13). 

(5) It is well known that the anisotropic factor A = pt—pr should vanish at the 
origin. However, in the present model the energy density under consideration 
has an inhnite central density. According to Misner and Zapolsky [36] one may 
use it as an asymptotic solution with a cut-off density below which the Eq. 
(27) will not valid. 

(6) In the present model, the stability of the matter distributions comprising 
of the anisotropic star has been attained. We have come to conclude this by 
analyzing the TOV equation which describes the equilibrium condition for 
matter distribution subject to the gravitational force (Fg), hydrostatic force 
{Fh) and another force (Fa) due to anisotropic pressure (Fig. 8). 

(7) For an overall view of the present study and also to have a physically viable 
model we put here the data which are available on some compact stars, e.g. 
Her X — 1 and SAX J 1808.4 — 3658 and Strange Star — AU 1820 — 30 (see 
Table 1 of the Ref. [IB])- The compactness factor u for these stars are 0.168, 
0.299 and 0.332 respectively. In comparison to these data we can conclude that 
our model represents a star with moderate compactness and also indicates 
that the star might be a compact star of Strange/Quark type rather than 
Neutron star. This is because, in general, superdense stars with mass-size 
ratio exceeding 0.3 are expected to be composed of strange matter [49] . 

Now according to Buchdahl [12], the maximum allowable compactness for a 
fluid sphere is given by ^ < |. Our compactness values for the chosen model 
(see Table 1) within this acceptable range and hence provides a stable stellar 
conhguration. 
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